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ABSTRACT 


The  concept  of  driving  force  per  unit  crack  front  or 
the  “force  tendency11  is  discussed*  This  is  a  more  general 
and  more  readily  useful  concept  than  the  Griffith  theory 
in  that  it  does  not  require  the  assumption  of  the  absence 
of  plasticity  and  there  is  no  complication  involving  fixed 
grips*  It  is  shown  that  the  force  tendency  11  exists 
for  any  stress  condition  and  that  it  can  be  determined  from 
strain  measurements  near  the  crack  head.  Various  formulae 
for  Jd  are  given  and  the  case  of  hydrogen  embrittlement  is 
discussed. 
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Onset  of  Fast  Crack  Propagation  in  High  Strength  Steel  and  Aluminum  Alloys 


G.  R.  IRWIN 

In  this  discussion  a  crack  extension^ force  tendency^  concept  will  "be  used. 
Since  applications  of  this  concept  to  fracturing  other  than  the  Griffith  theory 
are  not  widely  knovn  a  brief  description  of  it  will  be  given  first.  Measured 
values  of  the  crack  extension  force  at  onset  of  fast  fracturing  for  a  variety 
of  materials  and  test  conditions  will  be  shown  and  discussed.  Practical  use  of 
critical  values  of  crack  extension  force  will  be  illustrated.  The  paper  will 
conclude  with  a  few  remarks  on  hydrogen  embrittlement. 

The  force  tendency  is  defined  in  the  following  my.  Assume  one  has  under 
consideration  a  solid  object  acted  upon  by  external  forces  in  such  a  -way  that 
it  contains  a  growing  crack.  If ,  during  a  time,  gt,  the  work  input  from  these 
external  forces  is  S  and  the  change  in  stored  recoverable  strain  energy  is 
gU  we  can  write 

$ wE  (1) 

where  the  difference,  g  Wq,  must  consist  primarily  of  a  conversion  of  energy  to 
heat  within  the  object. 


Putting 


wg ,  .  jjl 

SWE=  £,  ?±  4  sJ±  ) 


(2) 


where  5-4  5-^®^  are  the  non-recoverable  and  recoverable  increments  of 


motion  of  Fi;  one  may  observe  that 

swQ  =  Z.  n  *  Zi  Fi 

$£±  -  SU 

(3) 

We  define  the  crack  extension  force  tendency  as 

*  .  - 

6v 

A 


in  the  limiting  case  as  S  A,  the  element  of  new  fracture  area,  approaches  zero . 
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The  only  reason  for  this  expression  to  differ  from  zero  i3  the  configurational 
change  of  the  object  due  to  extension  of  the  crack.  Its  value  depends  only  upon 
the  applied  forces,  the  configuration,  and  elastic  constants  of  the  material. 

Calculation  of  jb  for  the  central  crack  in  a  plate  subjected  to  simple 
tension,  cr  ,  shows  that  the  crack  extension  force  is  proportional  to  the  square 
of  times  the  length  of  the  crack,  as  one  would  anticipate  from  the  analysis 
by  A.  A.  Griffith (l) .  -Using  any  two  dimensional  situation  in  which  stresses 
are  applied  to  a  plate  containing  a  crack  along  the  x  axis  and  such  that  £xy 
along  the  x  axis  is  zero  (see,  for  example,  the  solutions  given  by  Westergaard) (2) 
one  can  readily  obtain  a  local  stress  interpretation  of  /&  .  If  the  Cartesian 
position  coordinates  of  the  stress  functions  are  replaced  by  polar  coordinates, 
r,  &  ,  using  one  end  of  the  crack  as  origin  of  coordinates,  as  r  approaches  zero, 

the  separational  stress,  <r~  ,  approaches  as  nearly  as  we  please  to 

«/ 

«rv  =  S22-4/Z  (1  +  sin  sin  2-S_)  (5) 

y  n  (2r)  *  2  2 

where  E  is  Young 5 s  modulus. 

Cracks  extend  by  a  process  of  development  and  joining  up  of  new  fracture 
origins  near  the  crack  tips.  In  connection  with  the  comparative  extension 
possibilities  of  cracks  of  different  length,  and  in  different  fields  of  overall 
nominal  stress,  or ,  the  result  of  the  preceding  paragraph  means  that  the 
development  of  new  fracture  origins  proceeds  under  the  influence  of  similar 
separational  stress  environments  for  two  different  cracks  when  the  values  of  ^ 
for  these  cracks  are  the  same.  In  general  for  cracks  of  various  length  and 
location  in  any  given  solid  material  the  value  of  Jb  measures  the  intensity  of 
the  crack  tip  stress  field  so  long  as  the  influence  of  plastic  deformations  ac¬ 
companying  fracture  extension  is  limited  to  the  close  neighborhood  of  the  crack. 


2 


i.(5 

The  omission  of  2f‘  F-  from  the  egression  for  the  force  tendency 

is  a  matter  of  choice.  For  example ,  when  we  analyze  the  longitudinal  yielding 
of  a  tensile  bar  with  no  crack  present  the  term  corresponding  to  is  of 
negligible  importance  compared  to  F  and  the  significant  component  of  the 

force  tendency  is  simply 


sF* 


z  F 


(6) 


On  a  per  unit  area  basis  this  is,  of  course,  just  the  longitudinal  stress. 
Fundamentally  ve  are  concerned  with  a  rate  problem.  ‘The  time  rate  of  con- 
version  to  thermal  energy  is,  approximately,  F  — 


Given  various 


values  of  the  force,  F,  one  wishes  to  know  the  corresponding  rate  of  yielding. 
This  must  be  found  experimentally.  The  familiar  result  is  that  a  plot  of 
£*  /  st  as  a  function  of  F  shows  the  yielding  rate  increases  from  slow  to 
fast  over  a  relatively  small  range  of  values  of  the  force  tendency.  (We 
neglect  work  hardening  effects  since  these  appear  to  have  no  counterpart  in 
fracturing.)  The  value  of  F  for  this  range,  per  unit  area,  is  the  yield 
strength  and  is  properly  regarded  as  a  significant  physical  strength 
characteristic  of  the  material. 

If  one  follows  an  analogous  procedure  in  the  study  of  fracturing,  various 
amounts  of  fracture  extension  forced  may  be  applied  to  a  test  specimen  con¬ 
taining  a  crack  and  the  response  measured  in  terms  of  time  rate  of  fracture 
extension.  It  is  found  generally  true  that  this  extension  rate  changes  from 
slow  to  fast  for  relatively  small  changes  of  the  crack  extension  force  •>  There¬ 
fore  the  critical  Jb  values,  thus  determined,  have  a  significance  relative  to 
resistance  to  fracturing  similar  to  that  of  yield  strength  as  a  measure  of 
resistance  to  yielding. 


i 

'  u 


may  be  calculated,  in  certain  crack 


Equations  by  means  of  which  ^ 
extension  situations  are  as  follows: 

1.  For  a  central  crack  in  a  larg^  plate  subjected  to  tension,  gr  , 
normal  to  the  crack  direction 

>  =  V  °~2a  (7) 

E 

where  &a  is  the  crack  length  and  E  is  Young's  Modulus.  This  is  the 
familiar  Griffith  expression.  The  same  expression  applies  to  an  edge 
crack  of  length,  a,  normal  to  the  tension,  <T  . 

2.  For  an  embedded  “disc  shaped"  or  “penny-shaped11  crack  in  a 
medium  in  which  the  stress  exerted  normal  to  the  plane  of  the  crack  is  cf  , 

*  -  h  -*4*-  (8) 

where  R  is  the  radius  of  the  crack  and  k^  is  k  (l  -72)/  7f  and  9  is 
Poisson's  ratio.  This  result  was  obtained  by  Sneddon(3)* 

3-  For  the  embedded  crack  as  above  but  with  the  addition  of  internal 

gas  pressure,  p,  Inside  the  crack 

Jb  -  kl  ( cr  4  p)£  R  (9) 

E 

Table  1  shows  values  of  ’*  ■J&jry  indicated  by  experiments  to  be  typical 
for  onset  of  fast  fracturing  in  sheets  and  plates.  There  are  undeniably 
large  influences  of  material  dimensions  as  indicated  by  the  values  for  24S-T4 
in  thin  sheet  as  compared  with  thick  plates.  It  is  also  interesting  to  note 
that  a  hard  and  somewhat  brittle  aluminum  alloy  75S-T6  is  tougher  than  mild 
steel  when  onset  of  fast  fracture  of  the  mild  steel  is  accompanied  by 
transition  to  predominately  cleavage  fracturing. 


Table  I 


Typical  values  of  the  crack  extension  force,  ■<*  , 
necessary  for  onset  of  unstable  fast  fracturing* 


Material 

Plate  of  Sheet 
Thickness 

Method 

Temperature 

Jj 

Steel,  ship  plate 
cleavage 

(in.) 

0.75 

a 

( Deg .  Cent . ) 

-  20 

(in.lbs/ln*1) 

100 

Aluminum  alloy  24ST4 

1.00 

b 

4  25 

4oo 

Aluminum  alloy  24ST4 

0.04 

a 

4  25 

600 

Aluminum  alloy  75ST6 

0.04 

a 

4  25 

300 

Aluminum  alloy  78ST6 

o.o4 

a  ! 

4  25 

150 

Polymethylmethacrylate , 
;as  cast  plates 

0.125  to 

0.50 

a 

4  25 

5; 

Polyesters,  plates 

0.125  to  0.50 

a 

4  25 

1 

Polymethyl-alpha- 

chloracrylate  plates 

0.125  to 

0.25 

a 

4  25 

2.8 

Vulcanized  natural 
rubber 

0.035 

c 

4  25 

15  to  30 

Glass,  lantern  slide 
covers,  moist 

0.02 

a 

4  25 

o.o4 

Glass,  lantern  slide 
in  2$  R  H  air 

0.02 

a 

4  25 

0.08 

Cellulose  acetate 

0.025 

a 

4  25 

14  ' 

Cellulose  acetate 

0.002 

a 

4  25 

28 

(a)  Centra!  crack  type  tension  specimens  in  which  the  unstable  crack  length  ms 
about  one  quarter  of  the  plate  width* 


(b)  Slow  bend  tests  of  notched  bars. 

(c)  Tear  tests  by  Rivlin,  Jnl.  of  Polymer  Scie.,  vgl.  10,  no.  3>  March  1953;  PP*  291-318. 
Various  compositions  give  values  of  Jb  from  10°  to  2  x  10*7  ergs/cm^. 


*Table  prepared  by  J,  A.  Kies  and  H.  L.  Smith 
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Mild  steel  and  plastics  such  as  Plexiglas  are  strain  rate  sensitive 
materials.  For  'uch  materials  considerable  variation  of  fa  can  occur 
dependent  upon  the  testing  procedure.  For  example  in  the  case  of  a  typical 
sheet  of  Plexiglas  for  a  given  fixed  magnitude  of  longitudinal  stress  the 
critical  length  of  central  crack  for  onset  of  fast  fracturing  can  have  three 
different  values  for  the  following  situations;  (a)  the  starting- crack  is 
represented  by  a  slot  made  with  a  jeweler's  hack-saw  prior  to  application  of 
the  load;  (b)  the  starting  crack  is  tapped  into  the  sheet  prior  to  load  ap¬ 
plication;  and  (c)  the  load  is  applied  first,  then  the  starting  crack  is 
formed  suddenly  by  an  impact.  The  values  of  fa  decrease  from  situation  (a) 
through  (c)  by  a  total  factor  of  about  2. 

On  the  other  hand  the  aluminum  alloys  are  much  less  strain  rate  sensi¬ 
tive  and  the  difference  in  fa  determined  by  methods  (b)  and  (c)  is  cor¬ 
respondingly  small,  G,  I.  Taylor  and  his  associates  shoved  that  the  dif¬ 
ference  in  static  and  dynamic  yield  strengths  of  steels  decreased  with  in¬ 
creasing  hardness.  Although  this  suggests  that  the  strain  rate  sensitivity 

h 

of  high  strength  steels  in  fracture  tests  may  not  be-, large,  some  reservations 
to  this  attitude  may  be  necessary  for  steels  which  are  temper  brittle  and  for 
embrittlement  by  hydrogen. 

R.  H.  Raring  and  J.  A.  Rinebolt  of  the  URL  Metallurgy  Division  have  been 
investigating  fracture  strength  of  high  strength  alloy  steels.  Their  ASM 
paper  at  the  1955  Metals  Congress  will  give  some  of  the  results  of  their  work. 
In  this  paper  use  will  be  made  of  several  of  their  test  results  and  of  test 
results  from  Syracuse  University  in  order  to  estimate  values  of  fa  for  cracks 
under  conditions  of  biaxial  constraint. 
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Figure  (l)  shows  a  fracture  obtained  with  a  tensile  bar  of  3AE  4340 
heat  treated  to  a  uensile  strength  of  260,000  p.s.i.  There  was  a  43$  reduction 
of  area.  The  region  of  the  origin  of  the  fracture  is  clearly  indicated. 

When  a  fracturing  process  changes  from  slow  to  fast  the  change  normally  oc¬ 
curs  rather  abruptly  with  the  more  slowly  fractured  region  differing  in  ap¬ 
pearance  from  the  region  cf  fast  fracture.  Although  IhiB  appearance  dif¬ 
ference  is  hardly  noticeable  in  materials  such  as  the  aluminum  alloys,  it  is 
a  common  feature  of  fracture  appearance  in  many  materials.  The  particular 
fracture  shown  on  Figure  (l)  shows  a  region  of  lighter  appearance  surrounding 
the  fracture  origin  and  ring-like  markings  near  the  periphery  of  this  region 
which  appear  to  be  hesitation  lines.  It  is  believed  the  position  of  change 
from  slow  to  fast  fracturing  in  this  instance  was  the  outer  edge  of  the 
zone  of  lighter  appearance.  The  expression  for  ^obtained  from  Sneddon's 
analysis  of  embedded  "penny- shaped"  cracks  can  be  used  to  calculate  the 
critical  value. 

In  this  example  the  average  longitudinal  stress  at  fracture  on  the  necked 
down  region  of  the  tensile  bar  was  3 50 >000  p.s.i.  and  the  measured  value  of  R 
was  0.025  in.  Thus 

=  120  in.  lbs/in2  (10) 

30  x  106 

Figure  (2)  shows  a  fracture  of  a  tensile  bar  of  the  same  steel  composition 
heat  treated  to  a  greater  hardness  than  would  be  customary  in  practical  ap¬ 
plications.  The  ultimate  tensile  strength  was  320,000  p.s.i,  and  the  average 
longitudinal  stress  on  the  slightly  necked  down  region  was  somewhat  less  than 
3-80,000  p.s.i.  If  one  assumes  the  lighter  toned  region  around  the  fracture 
origin  is  again  the  slow  fracture  zone,  calculation  as  above  gives 

-  150  to  200  in.lbs/in2 
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(ll) 


Fig.  1  -  Fracture  of  SAE  4340  Fig.  2  -  Fracture  of  SAE  4340 
tensile  bar,  U.T.S.  =  260,000  tensile  bar,  U.T.  S.  =  320,  000  psi 

psi  (see  first  example  of  Table  II)  (see  second  example  of  Table  II) 


Fig.  3  -  Fracture  of  SAE  4340  tensile 
bar,  U.T.S.  =  300,  000  psi  (see  third 
example  of  Table  II) 
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depending  both  upon  judgment  of  the  value"  of  R  and  upon  the  uncertain 
magnitude  of  the  load  drop  before  fracture  occurred. 

Figure  (3)  shows  another  tensile  bar  fracture  of  the  same  steel  again 
at  a  high  hardness  level.  The  ultimate  tensile  strength  was  expected  to  be 
about  300)000 'p.s.i.  Fracture  occurred  at  a  true  stress  of  293)000  p.s.i. 
at  a  surface  flaw  after  about  2$  reduction  of  area  had  occurred.  The 
fracture  appearance  shows  that  the  tendency  of  the  material  to  develop 
localized  plastic  shearing  near  the  bar  circumference  halted  the  spread  of 
the  flat  central  tensile  fracture.  Separations  and  weakening  incidental  to 
localized  plastic  shear  then  induced  a  change  in  the  direction  of  the  fracture 
onto  surfaces  of  maximum  shear  stress.  A  value  of  Jf  ms  computed  using 
equation  (7)  and  a  value  of  a  taken  from  measurements  on  the  photograph  to 
be  0.013  inches .  The  result  was  a  lvalue  of  117  in.  lbs/in2. 

The  results  of  Raring  and  Rinebolt  included  tests  in  which  sharply 
notched  tensile  bars  were  used.  A  larger  group  of  results  ms  also  available 
from  notched  bar. tests  at  Syracuse  University.  In  both  places  the  depth  of 
the  notch  was  chosen  so  that  the  holding  area  was  reduced  by  a  factor  of  two. 

A  fair  approximation  of  $  for  inward  extension  of  a  circumferential  crack 
can  be  made  using  the  following  expression 

Jf  =  -TTgi.2.3-  { 1  4  _  (  1  -  JL.  )l  (12) 

E  t  R0~a  ft  J 

Here  the  notch  depth  plus  the  average  depth  of  the  starting  crack  at  in¬ 
stability  is  a,  Rq  is  the  radius  of  the  bar,  and  O"  is  the  nominal  stress 
in  regions  of  the  bar  remote  from  the  crack. 

Using  the  work  at  Syracuse  University,  examination  of  the  notched  bar 
fractures  for  Various  size  bars  of  an  SAE  4340  steel  and  for  various  size  bars 
of  75S-T6  aluminum  alloy  was  made  with  the  purpose  of  selecting  a  situation 
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root  all  around  the  specimen.  The  graph  of  results  for  notch  strength  as  a 
function  of  size  shoved  in  each  case  a  plateau  region  for  the  small  sizes 
followed  by  a  drop  off  in  notch  strength  with  increasing  test  bar  diameter. 

In  both  materials  the  fractures  most  nearly  approximating  to  the  desired 
"inward  spread  of  the  crack"  situation  appeared  on  test  bars  at  the  beginning 
of  the  notch  strength  drop-off  region  of  bar  diameter*.  Thus  calculations 
were  made  for  a  0-35  inch  steel  bar  and  for  a  0.8  inch  aluminum  alloy  bar. 

For  the  SAE  4340  steel  bar,  using  the  value  of  notch  depth  for  "a"  in 
the  equation,  the  A  value  obtained  was  180  in.lbs/in2.  For  the  0.8  inch 
75S-T6  test  bar  the  calculation  on  the  same  basis  gave  100  in.lbs/in2. 

Table  II  gives  a  summary  of  critical  crack  extension  force  estimates  for 
the  examples  discussed  above.  It  will  be  noted  that  the  value  for  75S-T6  in 
Table  II  is  one  third  of  the  value  shown  for  this  material  in  Table  I.  The 
difference  is  due  to  the  additional  work  which  enters  into  the  fracturing 
process  for  a  crack  running  across  a  plate  in  regions  of  the  crack  close  to 
free  surfaces  of  the  plate.  Although  no  values  for  pertaining  to  running 
of  cracks  in  plates  of  high  strength  steel  have  been  measured  it  may  be  expected 
that  these  values  when  obtained  will  be  greater  than  200  in.lbs/in  .,  probably 
in  the  neighborhood  of  500  in.lbs/in2.  In  all  the  A  results  for  metals  pre¬ 
sented  above,  the  plane  of  the  fracture  is  perpendicular  to  the  direction  of 
greatest  rolling  or  forging  extension.  If  the  expected  weakest  fracture 
orientations  had  been  used  the  values  of  Jb  would  have  been  greatly  reduced. 

As  an  example  of  use  of  measurements  of  critical  Jr  values  consider  the 
case  of  fracture  propagation  which  appears  to  approximate  tne  unfortunate  Comet 
planes.  The  material  used  here  is  based  upon  work  by  Kies,  Smith  and  Brossman(4) 
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Table  II 


Critical  Crack  Extension  Forces  for 
Cracks  under  Biaxial  Constraint 


i 


n 


at  the  Naval  Research  Laboratory.  If  a  pressurized  fuselage  of  75S-T6  aluminum 


alloy  is  made  to  bear  service  hoop  stress  loads  of  20,000  p.s.i.  a  simple  cal¬ 
culation  using  300  in.lbs/in^  as  the  critical  lvalue  shows  that  a  crack 
longer  than  5  inches  will  be  unstable.  Now  it  may  be  shown  analytically  and 
ha3  also  be  demonstrated  in  NRL  laboratory  tests  that  a  crack  starting  from 
the  edge  of  an  opening  in  a  sheet  in  tension  is  subjected  to  a  crack  extension 
force  approximately  equivalent  to  that  of  a  crack  equal  in  length  to  the  full 
size  of  the  opening.  Thus  a  crack  started  by  fatigue  or  otherwise  at  a  window 
opening  in  such  a  pressurized  fuselage  would  propagate  unless  the  window 
dimensions  were  less  than  5  inches.  On  the  other  hand  a  more  cautious  practice 

in  plane  design  might  result  in  a  fuselage  of  24S-T4  aluminum  alloy  carrying 

2 

hoop  tension  stresses  no  larger  than  15,000  p.s.i.  Using  600  in. lbs/in  for 
the  critical  value  of  crack  extension  force  one  finds  that  in  such  a  plane 
cracks  developing  out  from  a  window  opening  would  not  become  unstable  until 
the  combined  length  of  crack  plus  window  opening  exceeded  18  inches. 

Numerous  practical  utilizations  of  the  crack  extension  force  concept  can¬ 
not  be  expected  until  supported  by  a  more  extensive  set  of  experimental 
measurements  than  is  currently  available.  The  limited  number  of  applications 
which  have  so  far  developed  pertain  to  the  "crack  running  in  a  plate"  type  of 
situation. 

Values  were  given  in  preceding  discussion  for  onset  of  fast  fracture  in 
situations  where  plastic  distortions  accompanying  crack  extension  are  re¬ 
stricted  as  for  an  embedded  crack  or  for  a  crack  at  the  root  of  a  sharp  notch. 
These  critical  crack  extension  force  values  were  rough  approximations  only. 

When  additional  experiments  provide  values  that  can  be  used  with  greater 
confidence  it  will  be  possible  to  calculate  the  depth  of  surface  cracks  and 
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diameter  of 'Eternal  cracks  which,  at  given  stress  levels,  may  be  expected  to 
cause  onset  of  fast  fractures . 

We  turn  next  to  the  influence  of  hydrogen  molecular  gas  pressure  within 
an  embedded  crack.  Carney,  Chipman,  and  Grant(5)  calculated  that  pressures 
ranging  from  105,000  to  220,000  p.s.i.  would  develop  at  room  temperature  in 
cavities  within  steels  containing  0.0001  to  0.001  per  cent  by  weight  of 

A 

hydrogen.  The  lower  figure  is  not  an  uncommon  hydrogen  content.  Selected 
regions  of  hydrogen  embrittled  steel  might  contain  0.001  per  cent  hydrogen. 

From  equation  (9)  it  would  appear  that  an.  internal  crack  capable  of  extension 
by  300,000  p.s.i.  tensile  load  under  conditions  such  that  hydrogen  pressure 
has  no  effect^  would  extend  equally  fast  at  200,000  p.s.i.  when  assisted  by 
100,000  p.s.i.  of  hydrogen  pressure.  However,  the  extension  rate  must  be  very 
slow,  allowing  time  for  diffusion  of  hydrogen  to  maintain  the  high  internal 
pressure,  in  order  for  hydrogen  to  have  its  maximum  influence. 

A  multiplicity  of  well  developed  cracks  only  one  or  two  of  which  are 

-7 

j 

Included  in  the  final  fracture  surface  is  a  typical  appearance  feature  of  i 

hydrogen  embrittled  steel  under  stress  rupture  test  conditions.  To  explain 
this  appearance  feature  one  needs  a  self -arresting  crack  extension  process. 
Keferring  to  equation  (9)  again,  one  may  note  that,  when  a  quick  increment  of 
crack  extension  occurs,  the  value  of  $  would  be  increased  by  the  increase  of 
crack  radius,  "a",  but  would  be  decreased  by  decrease  of  the  internal  pressure,  p. 
Whether  Jb  decreases  as  the  crack  extends  quickly  depends  upon  the  ratio,  p/cn, 
and  upon  the  initial  non-elastic  volume  of  the  crack.  Calculations  based  upon 
these  considerations  lead  one  to  conclude  that  minimum  internal  pressures 
ranging  from  £  (for  negligible  non- recoverable  crack  volume )  to  (5/10  cr‘ 

(for  relatively  large  non- recoverable  crack  volume)  are  required  for  the  crack 
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>'• 


extension  to  "be  self -arresting •  It  appears  therefore  that  the  mechanical 
internal  pressure  influence  of  hydrogen  in  steel  provides  a  plausible 
explanation  of  the  extensive  growth  of  cracks  as  a  self -arresting  process 
as  well  as  of  other  major  features  which  characterise  hydrogen  embrittlement. 
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